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Abstract 

We derive the energy of the surface between the normal and superfluid components of a mixed 
phase of a system composed of two particle species with different densities. The surface energy 
is obtained by the integration of the free energy density in the interface between the two phases. 
We show that the mixed phase remains as the favored ground state over the gapless phase in weak 
coupling. We find that the surface energy effects emerge only at strong coupling. 
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The issue of superfluidity of fermionic atoms in cold atomic (laser) traps has motivated a 
wide community as a consequence of the improvement of experimental techniques for dealing 
with these systems in the laboratory [l|, 0, H, 0, @, 0]- A richer physics emerges when 
the populations of the two candidates to form pairing are unequal. There will be unpaired 
fermions in the system. In this fascinating scenario, various alternatives for the ground state 
of these asymmetrical fermionic systems have been proposed in the literature [8|, 
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In Refs. [Tl|, 03 the energies of the mixed phase (MP) and the breached-pair superfluid 
phase [Io[ were compared for different density asymmetries of a fermionic gas composed of 
two particle species. It was shown that phase separation is the energetically favored ground 
state. This figure appears to gain strong experimental support [l5|, 13 . 

In these calculations the energy of the surface between the two equilibrium phases of the 
MP has been neglected, since this term is negligible in the thermodynamic limit (infinite 
volume), as considered in [ll|, 13] • 111 principle, the cost associated with the interface sepa- 
rating the phases could play an important role in finite atomic systems. For instance, the 
surface tension between the normal and superfluid phases can result in deformations of the 
minority component that are quit similar to those observed by the Rice experiments (l7j . 
However, the very poor knowledge of the surface energy function properties precludes a fully 
quantitative description of the phase separation stability in asymmetrical fermionic systems. 

In this paper we derive the surface energy that could be present in ultracold asym- 
metrical fermion superfluids, in order to verify the stability of the inho mog eneous phase 



separation [ll|, [12J over the homogeneous gapless superfluid (HGS) phase [10|, |13 



1. The Mixed Phase Total Energy 



We consider a dilute fermionic gas made up of a and b particle species (that could rep- 
resent, for instance, two hyperfine-spin states of fermionic atoms, such as 6 Li or 40 K) with 
densities n a and rib, m a "box" of volume Vbox- The MP state is formed by h a and n b (un- 



paired) atoms in the normal phase, immersed in a sea of 



BCS 
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BCS 



n (paired) atoms 



in the superfluid phase. Thus, the number densities n a and rib of the two atom species are 
accommodated within the box as n a ^ = xh a ^ + (1 — x)n. The parameters x and n are found 
by minimization of the energy of the mixed (inhomogeneous) phase ll|, [12J which will be 
defined below. 

For practical calculations we consider that the normal asymmetric state is in a bubble 
of radius R and the symmetric sup erfluid state is in the rest of the volume, defined by a 
Wigner-Seitz cell of radius Rw jl8| |. Then the total energy of the MP is given by: 



E MIX (n a , n b , a) 

V bu E N (h a ,h b ) + (V box - H«)E BCS l 



n) 



box 



xE N (n a ,n b ) + (l-x)E 



BCS, 



n 



(1) 
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where Vb u = ^wR 3 , x = V bu /V box is the volume fraction, a is the surface energy density, to 



be calculated below, E (n a ,n b ) 



(6tt 2 )5 
20tt 2 



-5/3 ^5/3 
"a I n b 



is the energy of the normal (unpaired) 
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particles, and E BCS (n) = ^ot^m ^(6tt 2 ^) 1 ^ 3 Aq is the energy of the superfluid state, 

where Aq is the mean field BCS gap parameter in the weak coupling approximation: 



,BCS 



Aa 

A = -^—e^, (2) 



e 2 



where /x BCS = /i BCS + /i BCS , kp = ^/2M/i BCS is the Fermi surface of the asymmetrical BCS 
superfluid, M = ™ a ™^ is the reduced mass, and a is the two body scattering length. As we 
pointed out earlier, for sufficiently large systems, meaning infinite volume, the last term in 
Eq. dl]) (oc R 2 ) is negligible in comparison with the first two ones (oc R 3 ). 

The limits found for the volume fraction are x = 1 if n b >> n a (the whole system is in 



the normal phase), and x = if n b = n a (the whole system is in the BCS phase) [ll|, [12 
Since in both limits the surface energy must be absent, we write the MP energy density 
E Mlx (n a ,n b ,a)/V box as 



E MIX (n a , n b , a) = E MIX (n a , n b ) + E s » (3) 

E MIX (n a ,n 6 ) 

= xE N (n aj n b ) + (l-x)E BCS (n) 

E su O) 

, A V 

+ [26(1/2 - x) + (1 - x)Q{x - 1/2)], 
R 

where Q(x) is the Heaviside "step- function" , defined as Q(x) = 1 if x > and Q(x) =0 if 
x < 0, E Mlx (n a ,n b ) is the energy density of the normal and superfluid phases and E su (cr) 
is the surface energy density. The lowest MP energy E MIX (n a , n b , a) is found by minimizing 
E Mlx (n a ,n b ) with respect to n and x, and the obtained x m i n is plugged in E su (<r) above. 



2. The Calculation of the Surface Energy 

We will now compute the surface energy between the normal and superfluid phases. The 
derivation is similar to the work of Reddy and Rupak [l8j which calculated, in the quark 
matter context, the surface tension between the normal and superconducting phases. To 
achieve this, they included the leading order gradient contribution to the free energy, whose 

(2) 

coefficient (termed there n 2 sc) was approximated as being independent of the chemical 
potential asymmetry and spatial variations of A. In Ref. [l8| determined by 

matching the mass of the fourth and eight gluon in the effective theory to the microscopic 
calculation while, here, we calculate the coefficients within the model under consideration. 

To describe the microscopic theory at the interface between the normal and superfluid 
phases, we begin writing the partition function: 



Z = J D[^ b ]D[^- s ^\ (4) 
where S(ip,ip^) = f dr f dx L>, with £ being the BCS Lagrangian 



3 



where ipi(x) describes a fermionic atom of the type i = a,b at x = (x,r), g > is the 
coupling constant, and /i aj b is the chemical potential of the non-interacting species. After 
equilibrium is reached the particles will have rearranged chemical potentials in each com- 
ponent of the MP. For this dilute system (5fi = flb ~^ a < /i afe ,and especially <5/i << A ), 
we consider that the local density approximation (LDA) holds, i.e. the chemical potentials 
in the individual phases of the MP are homogeneous and could be regarded as effective, 
including the contribution from the (external) trapping potential. We notice, however, that 
it appears that the LDA does not hold in the Rice experiments [l7| . and is consistent with 



the results of the MIT experiments |15l . 116 



We integrate out the fermions after introducing the usual Hubbard-Stratonovich field 
A(t, x), to obtain Z = J DADA*exp(— S e /f[A, A*]). The effective action is expressed by 

S eff [A,A*} = y!^£-TrlnG'- 1 [A(x)], (6) 



where J = ffdr J dx, and G 1 [A(x)] is the inverse of the Nambu propagator: 



G-^Aix) 



(7) 



A* (*) ar + £ + 

Near the critical temperature T c , where A/TKb can be considered a small quantity, the 



Ginzburg-Landau (GL) theory provides a successful description of the phase transition [19 
Although we are considering a zero temperature system, we have a physics motivated ap- 
proximation. Since the gap varies from Ao (on the superfluid side) to (on the normal side) 
at the interface between these two phases, we consider a slowly varying A and expand the 
action (in momentum-frequency space) up to forth order in A(q). At T=0 one should, in 
principle, include all terms of the expansion. However, if the spatial variations of A(x) occur 
on a length scale that is numerically large compared to Aq , it is a good approximation to 
keep only the leading term of the gradient in the expansion. This condition is guaranteed 
if V(A(x))/Ao << 1 [3], and we have verified that this is satisfied in the present work. 



This approach has been explicitly taken in Refs. [18|, |27j], where the solution found for A(x) 
justifies the approximation. 

S e// [A,A*] = 5>|A(g)| 2 + §|A(g)| 4 ], (8) 



I g ± L^in J (2tt 



d 3 k 
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where q = (q^iujm), with uj m = 2m7rT, A 

Q - n(<f, iuj m yj , u n = (2n + l)vrT, e^- + _. = e| = - fJ-b- As we are interested 

in the static limit of the corrections, we shall set iu m = 0. We also set m a = rrib = m. 
This brings about an enormous simplification to the problem without compromising the 
asymmetry of the system. The same consideration has been made in Ref. 21|. The coupling g 

d 3 k 1 



can be related to the (s-wave) scattering length a by ^ = J where = ^ |22| . 

Then we use this expression to regulate the ultraviolet divergence from the vacuum part of 
A. After performing the Matsubara frequency summation, the matter part of Il(g, 0) is 

,0, n f (ei )+n,(e b . ) 

written as IT m (cf,0) = / 44 k + m A k ~ m , where n f (x) = l/(exp(x/T) + 1) is the 

k + q/2 k-q/2 



4 



Fermi distribution function. Expanding n /( e g^/ 2 ) f° r \o\ << ^f' 6 ' with k F ,b = ^Jlm\i a f„ 
and taking the zero temperature limit we obtain 

f d 3 k © a - &(<1 2 /^-kq cos 9) + (a 6) 

n m fo o) = / — — F - - - - (9) 

W ' J 7 (2tt) 3 l/m[A; 2 + g 2 /4-m/i] ' 1 J 

where 6 a = Q(k F — \k\), 5 a = S(\k\ — k F ) is the delta function, and /i = \i a + The 
evaluation of the integral in the equation above yields .A(|<f|) = a + r y\q\ 2 + [J^i-J, where 



F 

b 

F^ n F 



a = N o (ifcTR - l): and 7 = Here we have defined, N = and k F = ^ 

is the "average" Fermi surface. This will allow us to study which phase is the most likely 
ground state as a function of the non dimensional interaction parameter 1//cf|o|. 

S S 

We now obtain the equation SA ^l = 0, and then take its Fourier transform to (x, 0): 



[a + /5|A(x)| 2 - 7 V 2 ]A(x) =0. (10) 

The solutions of the equation above depend on a, and can be divided into the limits a < 
and a > 0. 

A. n > 0, a < 

In absence of gradients, we obtain a space-independent solution that corresponds to the gap 
inside the superfluid region: 

|A|' = -f, (11) 

which has meaning only if the ratio % is negative. The trivial solution refers, obviously, to 
the gap inside the normal region. 

Seeking now space-dependent (one-dimensional) solutions, let us define a real function 

A(x) = Af{x) j20|, where A — \ \^\\ ■ Then Eq. ( TTUj) can be written for a < (as will be 
clear in the expression for the energy below, a stable theory needs f3 > 0) as 
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e^f(x) + f(x)-f(x) 3 = 0, (12) 

where £ 2 = A- . We find the solution of Eq. (TT2"]) 1 writing f (x) = 1 — e(x) , where e(x) << 1 
and keeping only linear terms in e(x): 

e( x ) = e-^ x/ t. (13) 

With the boundary conditions (B.C.) f(x — > oo) = 1 and f(x — > — oo) = 0, the satisfactory 
solution is f(x) = (1 — e(x))Q(x). Since inside the superfluid region A does not vary, A has 
to be identified with the BCS gap parameter. Namely, A = A . Thus the order parameter 
A(x) takes the form 

A(x) = A [l - e-^ x/ ^}Q(x). (14) 



Eq. (fT2|) also admits the solution f{x) = tanh(a; / 'y/2£)6(x), as well as the ansatz f(x) = i[tanh(:r/v2£) 
1], both satisfying the B.C. and exhibiting approximately the same profile at the interface. 
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Interface 
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x 



FIG. 1: Surface region between normal and superfluid components of the mixed phase. 



In Fig. ([I]) the gap profile is sketched as a function of x at the interface between the normal 
and superfluid phases. The energy density associated with the interface a la GL is 

E(A(x)) = aA(x) 2 + ^A{xf - 7 A(x)V 2 A(x). (15) 

From the equation above we see that deep inside the superfluid phase we have E(A(x — > 
oo)) = |Aq = Eo, and inside the normal phase, E(A(x —> — oo)) = 0. Since a stable 
interface requires (in order to fulfill mechanical equilibrium) equal energy densities on both 
the normal and superfluid sides, we must have that the energy inside the normal phase is 
also E Q . This is accomplished if we define a total energy density 



E T (A(x)) = E(A(x)) + e(-x)E , (16) 

The equation above clearly satisfies E T (x — > oo) = E T (x — > — oo) = E . Then, the surface 
(excess) energy density is defined to be: 

/oo 
[E T {A{x))-E ]dx, (17) 
■oo 

where the term subtracted in the equation above is the energy of the fully superfluid phase 
(or the fully normal state). If we multiply Eq. ffTUj) by A, integrate with respect to dx and 
substitute the resulting equation into the equation above, it is then easy to show that the 
(stationary) surface energy is given by: 

a = ^1 ( (°° [e(x) 4 + 0(-x) - l]dx\ . (18) 



2 

In Eq. (|18p we have eliminated (3 in favor of a. The x (numerical) integration yields 
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a 



251 r a A 2 , s 

24 ^-r- (19) 



Before pursuing further, it is worth mentioning that the surface tension could be obtained 
by pure phenomenological considerations. It is given by the product of the typical energy 
density difference between the superffuid and normal phases and the correlation length, as 
we have obtained. Nevertheless, the microscopic derivation serves to clarify and validate the 
results and allows direct extension to other systems. The result in Eq. (1191) is for a < 0, 
which implies positive surface energy and could, a priori, favor the HGS phase relative to 
the MP state. However, a < (or > — -) characterizes a system in the strong coupling 
limit which is beyond the validity of the mean field BCS approximation. In this range at 
least the leading order kpa corrections to the pressure and chemical potential should be 
included 2l[. We define: 



a + = * (l/kpa > -4/tt) = ^^M£A 2 , (20) 

where the true range of strong coupling in the equation above, such that the BCS part of 
the MP is a superfluid, is < l/k F a < 0. As we will see below, couplings beyond this 

point (i.e., l/k F a > 0), correspond to weakly coupled molecules in BEC 23|. 

In the case where a > or, in the weak coupling limit l/k F a < —4/ir (i.e., a small and 
negative), one would have to go to sixth order in the expansion in A. Then the energy 
expression, Eq. ( fl~5l) . is now 

E(A(x, 5)) = aA(x) 2 + ^A(x) 4 + -A(a;) 6 - 7 A(a:)V 2 A(x). (21) 

2 3 

The energy will have a stable configuration i.e., bounded from below, (for a > 0) provided 
the signs of the expansion coefficients are (3 < and 5 > 0, which in turn depend on 
the chemical potentials asymmetry. The space independent non-trivial equation for the 
minimum of Eq. (1211) is a — \/3\ Aq + SAq = 0, whose solution is 



A? 



l 0,± 



28 



1 ± v 7 ! - 4a5//? 2 l • (22) 



A real and positive solution also requires (3 2 > 4a5. In the point where a = 0, we have Aq + = 
0, corresponding to the normal solution, and Aq _ = \/3\/S representing a superfluid solution 

Thus the solution which gives the minimum for a > is Aq _ = j|j 1 + \J\ — 4a5/(3 2 
Proceeding as before, the surface energy is written as 



a a>0 



/oo 
[E T (A(x))^° - E^°] dx, (23) 
-oo 



To derive the expression above we have used the equation of motion [a + (3A(x) 2 + SA(x) 4 — 
7V 2 ] A (x) = 0, which yields 

E T (A(x))^° = |lA(x) 4 - ^A(x) 6 + e(-x)E^°, (24) 
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where Eq-° = ^aA^__ — ^Aq_ = ^Aq_ ^1 — Aq _J , due to the minimum condition. 
If we define a real function as A(x) = A 0t _g(x), the equation for A(x) is now 



cig\x) 



c 2 g(x)° + c 3 g{x)° - c 4 —g{x) = 0, 



' dx 2 



(25) 



where the coefficients c\ to C4 are all positive, and given by 



Cl 
Cl 



o:; 



(26) 



l/?|A^_; 
7- 



We write g(x) = 1 — ^(x), with x(x) << 1, and get after keeping only linear terms in x( x 

d 2 



dx 2 



X (x) = C lX {x) + C 2 



(27) 



where C\ = 5c3+ ^ 3c2 and C 2 = C2 ^ Ci ■ The equation above can be solved as x( x ) = 

y(x) — C2/C1, which yields y(x) = e~^^ x . However, C 2 = as a consequence of the 
minimum condition, and the final solution for the gap parameter satisfying the B.C. is 



A(x) = A _[l - e 



dx- 



Q(x) 



(28) 



The behavior of A(x)/A 0j _ is the same as the one shown in Fig{TJ Then, Eq. (1231) takes the 
form 



a 



a>0 



2 0> " 



^A(x) 4 - |a(x) 6 + e(-x)£^° - K-° 



POO 


g(*) 4 ( 


1 —00 





4 5 



(29) 



1 - -—A 2 



where we have used that Q(—x) 
above to obtain 



Q(x)dx, 

Q(x). We use the equation for A 2 , _ in the equation 



a 



a>0 



Wllere A = vfe 

< z < z 



9|/9| 3 e 



32 V^ 2 
9|/^| 2 V7 g(g) 



(s(u) 4 [l-p(^( M ) 2 ] -[l-^)]}du 



(30) 



where z m „ T < 1 will be find such that cr Q -° 



1 - e"", g(z) 



[1 + v 7 !^, with z = ±0-, 



0. Then we read the 
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limits for the interaction parameter from z in order to have a real surface energy: < 



a < Zmax/3 2 /45. In the second line of the equation above we have written y<5Ag_ — a = 

y/9(3 2 /85 g 2 (z)-p 2 z/45 = \0\/(2Vti) L(z), where L(z) = y/(9/2)g 2 (z) - z. Integrating in 
u we find that I{z) = -0.449 + l.eSSv 7 !^. Solving I(z) = we obtain z max = 0.924. For 
fixed <5/i, as is relevant to experiments, the parameters /?, 6 and 7 will also be fixed, and 
the only parameter of interest is 1/kpa encoded in a. Then we plot in Fig. [2] the function 

Ka a -°, where K = as a function of J0§tjzj l( z )i f° r z varying from to z max . The 
behavior of Ka a -° shows that the surface tension goes to zero for 1/kpa << — 1. 

B. jj, < 0, a > 

For completeness we analyze now the case a > in the BEC regime. We focus only on 
the space- independent solution of Eq. (fTOj) . After taking into account the normal energy, 
Eq. (fl~5|) at the minimum reads: 



E(A) = E N + -A 2 , (31) 



where E N = -±s&l(h - \fx\) 5/2 Q(h - (24|. In the BEC regime, where p. < 0, 



15 v^TT 



a 



V 2 in Ref. 0). In the case h < \p\, a « ^ (-^ - l). The point 



where a = gives t^— = 1. Substituting hp = kp+ J tF = V^rnii a +V 2 i!w± j n ^ s e q Ua tion, we 



k F a ^^^u^u^^g, 2 2 

find: 



v / /^+v / /^= 7S= • (32) 
V2m a 

Defining 

pi a = + /i, (33) 
Mb = I A* I - h, 
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we get for h < up to order h? 



E 



fa 2 



(34) 



where ^ = — t-. This is a well-known result for h = where, in this limit, tightly bound 
pairs with binding ener gy Eh, a nd nondegenerate fermions with a large and negative chemical 
potential are expected [22j, [24j. Eq. (134"I) serves, in this approximation, to define a critical 
value for the chemical potentials asymmetry, h c = = Then, < h < h c . 

The pair field A for fa < in this moderate coupling (1/kpa ~ 1) regime, where the 
mean field is not applicable, has to be determined by the number equation, 



n = — 



which gives, 



dE 



A 



A 2 da 



m 2 A 2 



n 2 



We plug in the equation above the total (free) density 

h a 3 h b 3 



h 3 

QlT 2 ' Q-K 2 67T 2 



n = n a + n b = 

where g(h/\»\) = (1 + Vl/"l) 3/2 + (1 - VM) 3/2 = 2 + | (ft)* + 0((VH) 4 ), to find 



2 



(35) 



(36) 



(37) 



A(fa<H) 



16 
37 



#(VH) I A* I 



(38) 



where \i is given by Eq. (I34p . From the equation above we see that the role of the chemical 
potential asymmetry is to lower the pairing field gap, as compared to the symmetric case 

(h = 0). Forh>Ua = £- a -^^ 



Lit 



yields 
E b 

The density now has the contribution of the normal energy: 

8E N A 2 da 2 m 3 / 2 



n 



dfi 2 0/x 3v^7r 2 



m 2 A 2 



Using n 



(2m^j,\) 1 ^ we g na jjy obtain 



A 2 (fa>H) 



8-^/2m|/x| 



3tt 



3/2 



2 m 3 / 2 



3^ 



87r-y/2m|/i| 



(39) 



(40) 



7T 



(41) 



where /i is given by Eq. ( |39i) . The chemical potential asymmetry for /i > also has the 
effect of lowering the pairing gap, as happens for h < \fj,\. The equation above sets the limit 
h cl <h< h c2 where h cl = fs and h c2 = (1 + 2 2 / 3 )f^. 
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The study of space dependent solutions of Eg. ({TO]) for a > 0, as well as the consideration 
of a possible surface energy between a phase with dilute free fermionic atoms, and a phase 



with condensed molecules within the BEC regime [24 , |25| , needs further investigation and 
will be considered in future work. 

In conclusion, we have obtained the surface energy that could arise in cold asymmetrical 
fermionic systems. The interaction limit, — A/tt < 1/kpa < (which implies a < 0, and 
consequently o > 0), shows that in the strong coupling regime (including the unitary limit 
i.e., infinite scattering length limit) a positive surface energy will emerge. We have shown 
that in the weak coupling BCS regimes, 1/kpa < — 4/V, where phase separation exists, 
it is found to be stable with no surface energy cost. In this regime, we found a lower 
limit for 1/kpa, which is simply an artifact of the truncation of expansion of the energy at 
sixth order. However, the vanishing of the surface energy for small kpa is an accurate and 
expected result. This is in complete accordance with the mean-field and Quantum Monte 
Carlo methods employed by Carlson and Reddy [2l|. Their results indicate that a HGS 
would exist only at very strong coupling, where A/fi > 1. They expect that in this regime, 
the surface energy cost could help stabilize the HGS phase over the MP state. 

We also found the expressions which fix /x in the BEC regime for h < and h > \fj,\, 
in terms of the molecular binding energy Ej,. We found that in both cases the role of the 
chemical potentials difference is to decrease the pairing gap size. 

We hope that the surface energy we have obtained could also be used in the theoretical 
investigations of the BEC-BCS crossover in asymmetrical fermionic systems (since in this 



entire regime of coupling, more coexisting exotic phases have been predicted to occur [26]), 
and also in the studies of superfluid properties of neutron matter. 



Note added: After completion of this work, a preprint appeared on arxiv [2j 
where there is investigated the role of the surface tension between the normal and superfluid 
regions of a trapped fermion gas at unitarity (kpa — > oo). They find in this regime that the 
surface tension is responsible for distortions in the atomic cloud in the boundary between 
these two regions. Although De Silva and Mueller use a different approach to obtain the 
domain wall energy at the superfluid-normal interface, they use the same expansion to 
calculate the coefficients a and 7. 
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